Free-carrier relaxation and lattice heating in photoexcited bismuth thin films 
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We report ultrafast surface pump and interface probe experiments on photoexcited carrier trans- 
port across single crystal bismuth films on sapphire. The film thickness is sufficient to separate 
carrier dynamics from lattice heating and strain, allowing us to investigate the time-scales of mo- 
mentum relaxation, heat transfer to the lattice and electron-hole recombination. The measured 
electron-hole (e — h) recombination time is 12-26 ps and ambipolar diffusivity is 18-40 cm 2 /s for 
carrier excitation up to ~ 10 cm' 3 . By comparing the heating of the front and back sides of the 
film, we put lower limits on the rate of heat transfer to the lattice, and by observing the decay of 
the plasma at the back of the film, we estimate the timescale of electron-hole recombination. We 
interpret each of these timescales within a common framework of electron-phonon scattering and 
find qualitative agreement between the various relaxation times observed. We find that the carrier 
density is not determined by the e — h plasma temperature after a few picoseconds. The diffusion 
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I. INTRODUCTION 

Intense ultrafast laser excitation of solids can alter dra- 
matically the interatomic forces by depositing substan- 
tial amounts of energy into free-carriers on a time scale 
that is short compared to energy relaxation. The result- 
ing non-thermal distribution of electrons in the bands 
of the material can drive coherent phonon motion, as 
well as dissipate its energy to the lattice by incoherent 
electron-phonon scattering. The evolution of the inter- 
atomic forces following photoexcitation is strongly de- 
pendent on the evolution of the electron distribution as 
it relaxes back to thermal equilibrium. 

Free carrier relaxation in metals, semimetals and semi- 
conductors is strongly influenced by the electronic den- 
sity of states at the energies where electrons and holes are 
present. Near equilibrium, electrons and holes in both 
metals and semimetals lie in an energy range of a few kT 
of the Fermi energy (k is the Boltzmann constant and T is 
the absolute temperature) . The significant difference be- 
tween the two is in the density of states at the Fermi level: 
for a typical metal, this is of the order of 1 eV _1 (per unit 
cell); for a semimetal like bismuth, it is of the order of 
10 -3 eV _1 . In metals, the relatively high density of final 
states gives rise to a momentum relaxation time of the 
order of a few femtoseconds 1 -, whereas the very low den- 
sity of states in semimetals suppresses relaxation towards 
thermodynamic equilibrium, particularly in equilibration 
between the electron and hole distributions 2 . Bismuth is 
an indirect semimetal — its electron and hole band ex- 
trema lie in different regions of the Brillouin zone (L and 
T points respectively) — so that electron-phonon scatter- 
ing is the dominant mechanism by which equilibrium is 



established between electrons and holes^. In semiconduc- 
tors, whereas thermalization of carriers within each band 
occurs on a sub-picosecond timescale, an energy band 
gap separates the electron and hole state energies, elimi- 
nating recombination by electron-phonon interaction and 
giving rise to typical recombination times of the order of 
ns. 

In the group V semimetals, the photoexcitation of 
modest plasma density leads to large amplitude atomic 
motion corresponding to a coherent zone-center fully 
symmetric (Ai 9 ) optical phonon 4 ' 5 . The initial pho- 
toexcited electronic distribution is determined by the 
intensity, photon energy and polarization of the inci- 
dent light. The intensity determines the overall num- 
ber of carriers, the photon energy their energy distribu- 
tion within the bands and the optical polarization gives 
rise to anisotropy in the distribution (i.e. differences in 
the photoexcited carrier occupation) between equivalent 
regions of the Brillouin zone. This distribution relaxes 
back to the thermal equilibrium distribution over time in 
a series of stages. Electron momentum relaxation (due 
to electron-electron and electron-phonon scattering) re- 
moves the anisotropy on a timescale of the momentum 
scattering time. The electronic energy distribution re- 
laxes more slowly because many scattering events are re- 
quired to change significantly the energy of an electron. 
When this process is due to electron-phonon scattering, 
heat is simultaneously transferred from the electronic to 
the vibrational degrees of freedom, raising the lattice 
temperature. Finally, scattering between conduction and 
valence bands, which can occur by electron-phonon scat- 
tering only when electrons are very near the bottom of 
the conduction band, establishes equilibrium in the occu- 
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pation of conduction and valence bands on the timescale 
of the electron-hole recombination time 2 -. 

Recently there has been considerable interest in the dy- 
namics of the Ai g mode in the limit of a dense electron- 
hole plasma, particularly for bismuth^—. Hase et al.— 
observed that the frequency of the mode is impulsively 
softened and chirped toward the equilibrium value over 
a few picoseconds. Murray et al.— demonstrated both 
experimentally using double pulse excitation and the- 
oretically using constrained density functional theory 
(DFT) and frozen phonon calculations, that the chirp 
was dominated by electronic softening and the evolution 
of the dense electron-hole plasma. Later, the theoretical 
studies were extended to include softening of the entire 
phonon spectrum^. The general results for the A\ g mode 
were confirmed by Fritz et al.— using femtosecond x- 
ray diffraction to measure the quasi-equilibrium position 
and curvature of the excited state potential. However, 
the results of measurement of ultrafast oscillations in the 
Debye- Waller factor by Johnson et al.— suggested that 
the coupling of the plasma with low-frequency acoustic 
phonon modes may not be well described in current the- 
oretical model a 17 ! 18 . 

The calculations in Ref. [ll| and [ID assume that a single 
electron- hole pair is created by each absorbed photon and 
that, whereas intraband scattering rapidly establishes a 
Fermi-Dirac distribution of carriers in each band, the 
electron-hole recombination time is substantially longer 
than the period of the relevant phonon and, as a result, 
the chemical potential for the conduction bands differs 
from that of the valence bands (two chemical potential 
model). This "two chemical potential" model is implicit 
in the previous analyses of low-temperature ultrasonic 
attenuation in semimetals^ (see also Ref. 0). In con- 
trast, Zijlstra et al.— performed frozen phonon calcula- 
tions on the zone-center mode, also assuming rapid scat- 
tering and equilibration throughout all conduction and 
valence bands, but using a single chemical potential, in 
which the carrier density is determined solely by the tem- 
perature of the electron-hole plasma. Johnson et al.— 
performed depth resolved femtosecond x-ray diffraction 
to measure the phonon dynamics and concluded that the 
two-chemical potential model may be appropriate only on 
a time-scale less than a single phonon period, transition- 
ing a single chemical potential model. Giret et al.— have 
analyzed the Bragg peak intensities, reported in Ref. H2I . 
using the model that assumes rapid thermalization of the 
carrier distribution throughout all bands and finds agree- 
ment with the experimental results that is as satisfactory 
as that obtained with the two chemical potential model. 
Recently Sciaini et al.— measured a 2-3 ps lattice ther- 
malization time of the photoexcited carriers in the limit 
of dense excitation using ultrafast electron diffraction, 
indicating a rapid energy transfer from carriers to the 
lattice. Still the dynamics of the photoexcited carriers 
are not well understood. The rapid transfer of energy to 
the lattice suggested in Ref. [2l| may call into question 
the detailed quantitative conclusions of analyses (includ- 



ing those of Refs. [IlTlllTlIllig) that have assumed 
insignificant heat transfer from the electron-hole plasma 
to the lattice vibration modes in the first few picoseconds 
following photo-excitation. 

Here we report results of ultrafast counter-propagating 
pump-probe measurements of the photoexcited carrier 
dynamics in a series of optically thick single crystal bis- 
muth films of various thickness at room temperature. By 
measuring the transient reflectivity signal on the oppo- 
site face of the film, we are able to separate the effects 
of carrier dynamics, lattice heating and strain, whereas 
the effects are coupled when probed on the same face 
as the pump2£. In this manner, we are able to mea- 
sure separately the ambipolar diffusivity along the trigo- 
nal axis and the carrier recombination rate. By compar- 
ing the reflectivity changes at the front and back of the 
film at delay times from 100 to 500 ps, we find that the 
film is heated much more at the front than at the back, 
consistent with recent grazing incidence x-ray diffrac- 
tion measurements on the thermal transport in thin Bi 
film a 23 i 24 and confirming unambigously that most of the 
heat transfer from the electron-hole plasma to the lattice 
occurs in the 5 ps following photoexcitation, before the 
plasma diffuses to a uniform density within the film2I. 

The principal aim of the current experiment is to esti- 
mate the timescales of momentum relaxation, heat trans- 
fer to the lattice, and electron-hole recombination. In 
particular, by investigating the diffusion of carriers to the 
back of films of various thickness, we estimate the mo- 
mentum scattering time. By comparing the final heating 
of the front and back sides of the film, we put lower lim- 
its on the rate of heat transfer to the lattice, and by 
observing the decay of the plasma at the back of the 
film, we estimate the timescale of electron-hole recom- 
bination. We interpret each of these timescales within 
a common framework of electron-phonon scattering and 
with reference to the known band structure of bismuth 
and find qualitative agreement between the various relax- 
ation times observed. We discuss the relationship of these 
results to early studies^ of low-temperature electron-hole 
recombination rates and to the relevant electronic den- 
sity of states in various regimes of electronic excitation. 
We argue, based on a picture of electron-phonon scatter- 
ing, that the relaxation of carrier momentum within the 
hole and electron bands of the initially excited plasma 
should occur at a rate typical of metals (with a relax- 
ation time of the order of a several femtoseconds), that 
cooling of the plasma and corresponding lattice heating 
should occur on a timescale of a few picoseconds and 
that electron-hole recombination should occur substan- 
tially more slowly than either of these processes. We 
argue that as the plasma cools towards room tempera- 
ture the electron-hole recombination rate and the plasma 
cooling rate should be of similar magnitude as a greater 
proportion of the excited carriers lie in the energy range 
where the conduction and valence bands overlap. 

The rest of this paper is organized as follows. In Sec- 
tion[H]we give the details of our experimental setup. The 
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results of experiments with low optical pump fluence are 
presented and discussed in Section Hill Measurements at 
high optical pump fluence are discussed separately in Sec- 
tion IIV1 In Section [V] we discuss the band structure of 
bismuth and the role of electron-phonon scattering in the 
various momentum, energy and electron-hole relaxation 
processes. In the concluding Section I VII we summarize 
the experimental findings and the principal conclusions 
of this study. Some technical details of electron-phonon 
scattering and carrier relaxation are presented in the Ap- 
pendix. 



II. EXPERIMENTAL DETAILS 

Various thickness single crystal c-axis Bi films were 
grown on 1 x 1 cm 2 sapphire substrates (c-axis ~ 0.5 
mm thick) by molecular beam epitaxy (MBE). We con- 
duct experiments using Tksapphire lasers; pulses of 85 fs 
duration, energy 3 nJ/pulse and wavelength centered at 
830 nm were produced at a repetition rate of 85 MHz 
(low density excitation) and pulses of 100 fs duration, en- 
ergy 4 fi J/ pulse, and wavelength centered at 800 nm were 
produced at a 250 kHz repetition rate (high density exci- 
tation). The laser pulses are short compared to the A\ g 
phonon frequency, 2.92 THz at room temperature. Ap- 
proximately 5 % of the pulse energy is used for the probe 
beam. A mechanical delay stage controls the pump-probe 
delay. The pump and probe beams are focused by lenses 
with 30 cm and 15 cm focal length respectively and have 
a diameter of ~1 mm before the lenses. The data are re- 
producible to better than 10% in the relative reflectivity, 
AR/R, to ^10~ 6 using a mechanical chopper, balanced 
detector and lock-in amplifier. 
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III. RESULTS AND DISCUSSION FOR LOW 
EXCITATION 



FIG. 1: (color online) (a) Ai 9 phonon oscillation. The inset 
shows the Fourier transform of Ai 9 oscillation, (b) Conven- 
tional surface-pump/surface- probe data. 



A. Surface pump and surface probe 



Fig. [T] shows ultrafast data for conventional surface 
(air/Bi) pump and surface probe for four different thick- 
ness films, ranging from 185-385 nm. We estimate the 
absorbed fluence to be about 0.8 /xJ/cm 2 . The reflectiv- 
ity data shows oscillations at 2.88 THz, near the room 
temperature A± g frequency, demonstrating that we are 
in the low fluence regime where the phonon frequency is 
not substantially softened by the photoexcitation. Phe- 
nomenologically, the relative reflectivity change is 



AR(t) dhiR dhiR 
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to first order in the A\ g phonon coordinate x, pho- 
toexcited carrier density n, strain 77, and lattice tem- 
perature T, all of which are functions of time t since 
laser-excitation. Generally the partial derivatives are not 



known; however, we can deduce their signs from the data 
as discussed below: dR/dx > (x driven towards the 
symmetric non-Peierls distorted structure), dR/dn > 0, 
dR/dr) < for tensile and > for compressive strain, 
and dR/dT < 0. We note that, in expanding the reflec- 
tivity as in Eq. ([T]), we do not presuppose any particular 
constraint on the carrier density n — for example, under 
certain physical conditions, n could be determined by the 
electronic temperature, T p . We assume second order and 
higher terms are negligible. 

The difference between the surface probe data (Fig. 
OJb)) for various film thicknesses is small. This is to 
be expected; because of strong absorption, carriers are 
generated near the surface (~ 15 nm laser penetration 
depth) and the probe beam only detects the reflectivity 
change in this region. We see AR(t)/R first increasing 
rapidly due to the sudden increase in carrier density and 
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FIG. 2: (color online) Counter-propagating surface- 

pump/interface-probe data (solid line) and fits for carrier sig- 
nals (dash line from 0-80 ps), described in the text. Each 
set of data is collected up to the time when the first acoustic 
strain is observed. 



then slowly decreasing as the carriers relax via multiple 
processes. As discussed above, the oscillations at early 
times are due to the A lg phonon ( clearly resolved in Fig. 
Ha) but aliased by the relatively coarse sampling in Fig. 
Hlb)). The overall reduction of the reflectivity signal dur- 
ing the first few picoseconds is due to a competition of 
effects from carrier relaxation through ambipolar diffu- 
sion, recombination, lattice heating and the propagation 
of an acoustic pulse generated through a combination of 
rapid thermal expansion and the acoustic deformation 
potential interaction. The signal becomes negative after 
approximately 10 ps as the strain pulse leaves behind a 
thermally expanded region at the surface, which decays 
on the much slower time scale of heat conduction in the 
lattice. We note that an optically thin sample would 
mitigate the effects of the inhomogeneous carrier distri- 
bution; however it would not decouple the effects of the 
coherent phonon, recombination, lattice heating or strain 
on the reflectivity change since they occur on overlapping 
time scales. 



B. Surface pump and interface probe 

Fig. [2] shows the relative reflectivity change for the 
probe incident on the Bi/sapphire interface. Clear differ- 
ences are seen for various film thickness, compared to the 
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FIG. 3: (color online) The peak time, t p , for carriers to arrive 
at the interface of the film as a function of film thickness, L, 
showing an approximately linear relationship with a slope on 
order 30 ps//im. 



surface-probe data. There is no coherent Ai g phonon sig- 
nal present on the interface-probe data, consistent with 
a localized excitation of the zero group velocity mode 
in the laser-excited volume of the optically thick films 25 . 
Thus, to the extent that we can ignore heating due to 
recombination at the interface, both the first and fourth 
terms on the right hand side of Eq. ([T]) are absent in 
the interface signal. The bipolar feature at late times 
(> 50 ps) is attributed to the acoustic strain pulse. We 
use the arrival time of the acoustic pulse to calibrate the 
film thicknesses relative to the thinnest film (185 nm) 
which we measured using grazing incident x-ray reflectiv- 
ity. The measured speed of sound for the thinnest film, 
2.15 ± 0.2 km/s, compares with 1.972 km/s of Ref. [H 
The uncertainty arises from the difficulty in determining 
the precise pulse propagation time due to the finite tem- 
poral width of the acoustic pulse. Well before the arrival 
of the acoustic pulse, the signal peaks at a time (t p ) that 
is approximately linearly proportional to the film thick- 
ness (see Fig. [3]). This peak is attributed to carriers that 
have diffused across the film before recombining. Note 
that the carriers' signal is well separated from the acous- 
tic pulse which takes approximately 15t p to arrive at the 
interface. 

The shape of the early-time component (carrier peak) 
of the reflectivity signal depends on both ambipolar dif- 
fusion and recombination. To model the dynamics, we 
assume that carrier diffusion is one-dimensional (in the 
coordinate z along the c-axis). This is appropriate since 
the laser pump diameter, ~150 /im, is many orders of 
magnitude greater than the film thickness. In the limit 
that the recombination time r and the ambipolar diffu- 
sivity D are independent of carrier density, n, 
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FIG. 4: (color online) Comparison of 185 nm film between 
surface-probe and interface-probe data both pumped from Bi 
surface at rig ~ 10 20 cm~ 3 assuming linear absorption. Mul- 
tiple acoustic pulses and echoes are shown. 
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D (cm 2 /s) 


t (ps) 


185 nm 


18±5 


12±2 


275 nm 


24±6 


26±4 


305 nm 


28±5 


24±3 


385 nm 


40±14 


26±6 



TABLE I: Fit results from various films. 



It is instructive to consider the limiting case of an in- 
finitely thick sample with a 5-function excitation at t = 0, 
for which 
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and we find that the maximum density at a given depth 
z occurs at time, 



(4) 



Note that deep enough into the bulk t p is linearly pro- 
portional to z and 





> y/Dr/A. 



(5) 



This approximation is reasonably accurate for the inter- 
face z — L of a relatively thick film. Thus, the linear 
dependence seen i n Fig. [3] suggests that our samples sat- 
isfy L 3> y/Dr/A. Moreover, Eq. ([5]) allows us to make 
a rough estimate for D ss 70 cm 2 /s given the ~30 ps//xm 
slope in Fig. [3] and ~25 ps decay from the reflectivity 
signal in the time period 10-50 ps after the pump. 
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FIG. 5: (color online) Simulation for 185 nm film using diffu- 
sion constant 18 cm 2 /s, recombination time 12 ps, and lattice 
thermalization time 1.5 ps for (a) carrier density (solid line 
with left axis) and energy density (dash line with right axis) 
and for (b) the increased lattice temperature as a function 
of depth at various time delays with initial carrier density 
n =2xl0 18 cm" 3 . 



To better model the data (and allow for variations in 
the individual samples), we solve Eq. ([2j numerically 
for the finite crystal case, treating D and t as fitting 
parameters. We assume that the carriers are confined at 
the surface (z = 0) and bismuth-sapphire interface (z = 
L), dn/dz\ z=0 L = 0, which is valid so long as surface 
recombination is slow. The initial photoexcited carrier 
distribution is assumed to decay exponentially over a 15 
nm laser penetration depth with a single electron-hole 
pair per incident photon. The initial carrier density no 
for the data in Fig. [5] is estimated to be about 2xl0 18 
cnr 3 at the surface. The probe is distributed across 
the laser penetration depth. The overall amplitude is 
normalized to the individual data of various films. The 
results of the fits are summarized in Table HI 

From the ambipolar diffusivity, D, we can make a 
rough estimate of the momentum relaxation in the re- 
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laxation time approximation. Assuming an average band 
effective mass m* of a free electron, and taking the relax- 
ation time as r m*D/(kT p ), we obtain an estimate of 
the scattering time, equal to 0.6 ps for a room tempera- 
ture plasma. In the initial hot photo-excited plasma, for 
which kTp w 0.5 eV, the relaxation time reduces to ap- 
proximately 30 fs. This is reasonably consistent with the 
results based on density of states arguments in Section 

El 

We note that the height of the first peak in AR/R 
for the low-fluence data indicates that the rate of carrier 
recombination prior to t p is not substantially faster than 
in the time after t p . If it were, a substantially smaller 
peak in AR/R would be observed. As we will see later, 
this is not the case for high-fluence data. 

In Fig. |4l the surface pump-probe signal drops be- 
low zero after about 10-20 ps, with a slow recovery on 
the nanosecond time-scale. This initial drop is caused by 
acoustic strain propagation from the surface through the 
optical absorption depth, leaving a thermally expanded 
lattice in its wake, while the subsequent very slow de- 
cay is attributed to cooling due to heat conduction. We 
note that the background reflectivity change due to this 
thermal expansion is much greater at the surface of the 
film than the corresponding background at the interface. 
Therefore the front surface of the film has been heated 
much more than the back interface. This is consistent 
with recent time-resolved x-ray diffraction experiments 
on the thermal transport in thin Bi, in which the thermal 
gradient is maintained for a time on order nanoseconds 
for comparable film thicknesse s 23 ' 24 . This demonstrates 
that lattice heating must occur substantially faster than 
the time (~ 5 ps) taken for diffusion to produce a uni- 
form plasma in the 185 nm thick film. Thus, the carriers 
cool by heat transfer to the lattice in a time substan- 
tially less than 5 ps and the carriers reaching the back 
surface of the film must be relatively cold, compared to 
the initial photo-excited plasma. This is also consistent 
with electron diffraction measurements of thin films re- 
ported in Ref. [2l| that suggest a lattice heating time of 
2 — 2.5 ps. Although the plasma cools in substantially 
less than 5 ps, its density decays with a much longer time 
constant (10-30 ps). 

This indicates that the decay of the plasma density is 
determined by the electron-hole recombination rate and 
we therefore estimate the recombination time to be of the 
order of 20 ps. (However, bearing in mind the discussion 
of cooling of the hot plasma, below, we cannot entirely 
rule out the possibility that the cooling rate of the plasma 
decreases as it cools towards the lattice temperature.) 
Using separate decay rates for carrier recombination and 
lattice thermalization: 

dn(z,t) _ n(z,t) | ^d 2 n 
dt r e - h dz 2 ' 

dE(t) _ E{t) 



= nwe -t/n I( J_ + L )n (z,t), (8) 

dt T e - h Tq 

where T e -h is recombination time, E(t = 0) = hut is 
photon energy absorbed per electron-hole pair with decay 
time tq after excitation, C v is lattice specific heat, and 
T; is lattice temperature, we show carrier density, energy 
density and lattice temperature change (without taking 
thermal diffusion into account) in Fig. [5] for the 185 
nm. The hot plasma (t — 0) cools much faster than the 
electrons and holes recombine. 

The rate of temperature loss in the semimetal plasma 
is comparable to the typical rates observed in metals (a 
few picosecond or less). However, the rate of electron- 
hole recombination (or, equivalently, scattering between 
the T and L valleys) is substantially slower. Our values of 
the recombination time, as observed by the decay of the 
plasma at the back side of the film during the interval 10- 
50 ps following photo-excitation, are considerably larger 
than a simple extrapolation of the fit of low-temperature 
acousto-magneto-optical effect data, given in Lopez—. 



IV. RESULTS AND DISCUSSION FOR HIGH 
EXCITATION 

For the carrier density > 10 20 cm -3 , the dynamics be- 
come more complex. Fig. [6] shows data as a function 
of time at a range of pump laser fluence for (a) surface 
pump and probe of 185 nm film and (b) surface pump 
interface probe of 385 nm film , scaled in peak amplitude 
for comparison. The individual inset shows how the am- 
plitude varies with fluence for (a) 185 nm film and (b) the 
different thickness films. The peak intensity of the carrier 
signal in the surface-pump/interface-probe data, shown 
in Fig. [6^b) , saturates and the peak time t p becomes 
longer and carrier recombination becomes faster. On the 
other hand, for the surface pump and probe, the max- 
imum (carrier density) and minimum (thermal) compo- 
nent of the relative reflectivity remain linear up to much 
higher fluence. We were not able to determine the level 
at which this signal saturates due to potential artifacts 
at the highest fluence for the surface probe. Nonethe- 
less, it appears that at high densities the recombination 
and diffusion become nonlinear. The increased plasma 
decay rate in the initial dense plasma reduces the am- 
plitude ARnfjnaxj/R of the carrier peak of the counter- 
propagating reflectivity at t — t p by up to a factor of 5 
at the highest fluences, indicating that the plasma den- 
sity decays at an average rate of 0.3 ps -1 at early times 
(t < 5 ps), comparable to the rate of heat transfer to 
the lattice. There are a few possibilities for the increased 
recombination in the non-linear regime: (a) the Auger 
channel begins to dominate over the phonon recombina- 
tion channel at high enough density (when n > y/^y/A, 
where 7 is the linear recombination rate and A is the coef- 
ficient for Auger recombination); (b) the phonon-assisted 
recombination rate is enhanced due to the presence of 
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FIG. 6: (color online) (a) Surface-pump/surface- probe data 
for 185 nm film at various absorbed fluence. The inset shows 
the maximum intensity at different fluence, which is similar 
for various film thickness, (b) Surface-pump/interface- probe 
data for 385 nm film at various absorbed fluence. no ~ 3x 
10 19 cm -3 is estimated from absorbed fluence of ~11 uj/cm 2 
assuming linear absorption. The peak time is delayed for 
higher fluence. The inset displays the peak intensity for var- 
ious film thickness at different fluence. From lower to upper 
lines are 385, 305, 275 and 185 nm respectively. 



more thermally excited phonons because the lattice is it- 
self heated; (c) the opening of new recombination chan- 
nels because the electronic structure is substantially al- 
tered when the Peierls distortion is partially reversed in 
the photo-excited material — for example, if the band 
overlap increased, this would open up a higher density of 
states for the indirect phonon-assisted recombination, or 
if the direct gaps either at T or L were reduced enough, 
we could have phonon-assisted recombination across the 
direct gap (see Section fVl below) . At this point, the na- 
ture of the nonlinearity cannot be determined without 
further experiments and more sophisticated analysis. 



V. DISCUSSION — BAND STRUCTURE AND 
ELECTRON-HOLE RECOMBINATION IN 
BISMUTH 

Bismuth has two atoms in the primitive cell and 5 va- 
lence electrons per atom, giving 10 valence electrons per 
primitive cell, which almost completely fill the lowest five 
valence bands. However, due to a 38 meV band overlap 
between the fifth and sixth bands, electrons occupy the 
sixth band near the three L-points in the Brillouin zone, 
with a carrier density ~ 10 17 cm~ 3 (~2x 10~ 5 per unit 
cell) at low temperature^ (see also Ref. H3). A corre- 
sponding number of holes occur near the T-point in the 
fifth band. The Fermi level lies approximately 26 meV 
from the conduction band minimum and 12 meV from 
the hole band maximum. The electronic density of states 
per unit cell is approximately 2 x 1CP 3 cV _1 for elec- 
trons and 3 x 10~ 3 eV" 1 for holes at the Fermi energy. 
This gives an estimate of the typical density of states in 
the region of the overlap of the conduction and valence 
bands. On the other hand, the typical density of states 
is approximately 1 eV _1 at energies 0.5-1.0 eV from the 
Fermi level^S. Fig. [7] shows schematically the electronic 
bands and scattering processes relevant for momentum 
relaxation, heat transfer to the lattice, and electron-hole 
recombination. 

Phonon assisted electron-hole recombination in Bi was 
discussed by Lopez in 1968. 3 From the acoustomagneto- 
electric effect, Lopez measured an electron- hole recombi- 
nation time ~10 -8 s below 4 K, which decreased rapidly 
above 6 K, assisted by the absorption or emission of 
phonons. The acoustic attenuation is fitted with a com- 
bination of (temperature independent) defect scattering 
and two activated processes, with characteristic temper- 
atures of 43 K and 130 K, which are consistent with 
recombination, assisted by two distinct phonons of en- 
ergy Tiuj equal to 3.7 meV and 11.2 meV, respectively, 
corresponding to acoustic and optical phonons at the X- 
point of the Brillouin zone. Lopez considered the assign- 
ment of the component with activation temperature of 
130 K to electron- hole recombination assisted by the op- 
tical phonon at X to be tentative and allowed for the 
possibility that this attenuation mechanism could be as- 
sociated with hole scattering between valence band ex- 
trema at T and elsewhere in the Brillouin zone. 

Extrapolation to room temperature of the expression 
given in Eq. 26 of Lopea^, gives an electron-hole recom- 
bination time of approximately 0.5 ps, dominated by the 
11.2 meV (130 K) phonon (if one neglects the component 
with activation temperature of 130 K, a room tempera- 
ture recombination time of 6 ps is obtained). However, 
although Lopez begins with a very general expression for 
carrier relaxation due to electron-phonon scattering in 
Section V.A of Ref. i, the further derivation leading to 
Eq. 26 and the detailed analysis of the experimental data 
depend on the assumption that the system is at very low 
temperature and that the density of electronic states does 
not vary significantly over energies of a few kT. Clearly 
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this assumption breaks down for Bi at room tempera- 
ture, where kT is comparable to the valence-conduction 
band overlap, and is entirely invalid for a photo-excited 
electronic plasma at temperatures, where kT ~ 0.5 eV. A 
quantitative prediction of the room temperature electron- 
hole recombination time from the results of Lopez is ques- 
tionable, but values of the order of 0.5-10 ps are not un- 
reasonable. 

In Appendix A, we give a more general discussion 
of the role of electron-phonon scattering in establishing 
equilibrium of the chemical potentials of various bands 
(including electron-hole recombination as a particular 
case) and in transferring heat from the electronic plasma 
to the lattice. The phonon-assisted rate of relaxation of 
the chemical potential difference between two bands, 1 
and 2, can be expressed qualitatively in terms of typical 
electron-phonon matrix elements and the average density 
of available final states, to which carriers can scatter by 
phonon absorption or emission: 



2tt 



71- 



{\M\ 2 )N(Tiu){D avall )/n 



(9) 



where (|-M| 2 ) is the average-square electron-phonon ma- 
trix element, which is expected to be of similar magnitude 
for short-wavelength acoustic modes and optical modes 
throughout the Brillouin zone, N{TwS) is the phonon oc- 
cupation number, (D ava u) is the density of final elec- 
tronic states per unit volume available at energies eihu, 
where e is averaged over the energy of carriers in the ini- 
tial band and u> is the typical optical phonon frequency, 
and n is the carrier density. The effective density of avail- 
able states includes a blocking factor (1 — / for electrons 
and / for holes, where / is the occupation of the final 
state) arising from the Pauli exclusion principle. Note 
that in our room temperature experiments all the bis- 
muth phonon modes are in the classical regime, with 
N(Hlo) ~ kTi/hco, where 7] is the lattice temperature. 
The total carrier scattering rate j p is a sum over the 
rates to all bands, including intraband scattering. We 
show in Appendix A that the rate of the heat transfer to 
the lattice from a high-temperature electron-hole plasma 
gives a plasma cooling rate, 



7Q 



ti'jj 



(10) 



where T p is the plasma temperature, 7) is the lattice 
(phonon) temperature and u is the typical frequency for 
phonons involved in the electron-phonon scattering pro- 
cess. We will discuss below the consequences of these 
results for electron-hole recombination and temperature 
relaxation in various regimes of photoexcitation. 



A. Momentum and energy relaxation of the hot 
plasma 

Immediately after photo-excitation, electrons and 
holes are created with energies approximately 0.75 eV 



into the conduction and valence bands, respectively. At 
these energies, the electronic density of states in each 
band is of the order of 1 eV -1 and, based on Eq. [9l we 
expect a momentum relaxation rate "f p in this hot plasma 
approximately 10 3 faster than the momentum relaxation 
rate in a room temperature plasma. Taking our crude 
estimate, above, for the carrier relaxation rate in the 
room temperature plasma, based on extrapolation of the 
low-temperature analysis of Lopez, where we expect a 
momentum relaxation time of the order of 0.5-10 ps, we 
would anticipate momentum relaxation time in the highly 
excited plasma to be of the order of several femtoseconds. 
Indeed, this is in good agreement with the momentum re- 
laxation time inferred in recent experiments comparing 
spontaneous and impulsive stimulated Raman scattering 
in Bi and Sb, corresponding to the decay rate of the elec- 
tronic forces driving the low-symmetry E g mode . In 
these experiments, a decay time, 1/7 P ~ 5 — 10 fs is ob- 
served for the Eg force at room temperature for Bi2£. On 
the other hand, taking a typical optical phonon energy 
in Bi to be approximately 10 meV, we can use Eq. ITOlto 
estimate a temperature decay time in the highly excited 
plasma to be 



Tiuj 



m 

TlLU 



1 _ 

7 P ~ 



4 ps. 



(11) 



This estimate is in good agreement with inhomoge- 
neous lattice heating observed here and with the elec- 
tron diffraction experiments reported in Ref. [2l] and is 
consistent with the diffusion of a cold plasma seen here. 

As the plasma cools, it is difficult to reliably esti- 
mate the change in the plasma cooling rate: Taking 
an estimate of the electron density of states based on a 
Dirac model of the bands near the L-point would imply 
7Q ~ an estimate of the density of states based on 

a parabolic band model (more appropriate for the holes 
near the T-point) would imply ~ T p In real- 
ity, we expect the temperature dependence of 7q to be 
relatively weak, with the changes in the factors (D ava u) 
and Tiui/kT p almost canceling each other as T p becomes 
smaller. 



B. Electron- hole recombination in the hot plasma 

For the hot plasma (with kT p ~ 0.5 eV), assuming 
phonon-assisted electron-hole recombination between the 
electrons near the L-point and holes near the T-point 
in the Brillouin zone, we expect two factors to strongly 
suppress the recombination process, compared with the 
intra-band scattering responsible for momentum relax- 
ation: (1) the small energy overlap of the conduction and 
valence bands implies that most electrons are at energies 
where they cannot scatter to the valence band (and sim- 
ilarly for holes scattering to the conduction band); (2) 
the density of states in the conduction-valence overlap 
region of the bands is much smaller than that typical of 
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FIG. 7: Schematic of the Bi band structure near the Fermi 
energy, showing typical electron-phonon scattering processes. 
There is a small, ~38 meV overlap between the valence band 
maxima at the T point (holes) and the conduction band min- 
ima at the L points (electrons), such that carrier recombina- 
tion occurs over a small region of phase space involving the 
emission or absorption of a phonon with momentum near X. 



the hot carriers within the band, which determines the 
momentum relaxation time. Including both factors (each 
approximately 10~ 3 ) suggests that the phonon-assisted 
electron- hole recombination is a factor of ~ 10~ 6 smaller 
than the momentum relaxation time. This rate is so 
small that this recombination mechanism is effectively 
closed off and it is likely that the actual electron-hole re- 
combination rate is determined by other processes, most 
likely Auger recombination, which would be expected al- 
ways to dominate at high plasma densities. 



C. Relaxation in the cool plasma 

Momentum relaxation, cooling and carrier recombina- 
tion in the plasma that has cooled, close to room tem- 
perature, is determined by two different aspects of the 
density of states. Only interband scattering between the 
electron and hole bands contributes to the electron-hole 
recombination rate, "f e -h, whereas all bands contribute 
to the momentum relaxation rate. Thus, it is clear that 
7e-ft < 7p- The momentum relaxation and the cool- 
ing rate of the lattice are still related by Eq. [TU1 but 
T p w Ti and the ratio Jq/j p is of the order of 1/10. It 
is more difficult to estimate the relative size of 7q and 
the electron- hole recombination rate, J e -h> although one 
would expect them to be of similar order of magnitude. 
To determine which is larger would require a more pre- 
cise knowledge of the electron-phonon matrix elements 
for intra-band and inter-band scattering than is currently 
known. 



VI. CONCLUSIONS 

We summarize the essential observations of the exper- 
iment: 

1. Most of the excitation energy of the electron- hole 
plasma is transferred to the lattice before the plasma 
becomes uniform in the films. We know this because 
the residual AR/R after 150 ps is much larger at the 
front face than at the back face. This residual dR/R 
can only come from heating of the lattice. This means 
that much more energy was transferred to the lattice at 
the front of the film than at the back. (Heat conduction 
through the lattice is too slow to significantly equilibrate 
the temperature in the film in 150 ps.) 

2. The plasma becomes uniform in approximately 5 ps. 
This comes from an analysis of the short time (<10 ps) 
back-side reflectivity signal, based on a generic picture 
of a plasma diffusing from the front-side of the film, as 
shown in Fig. [5] (a) using EqJS] This signal can only 
come from the electron-hole plasma since strain and heat 
conduction from the front of the film are too slow to reach 
the back side of the film in less than 50 ps. 

3. The back-side reflectivity AR/R decays with a de- 
cay time of approx 20 ps during the 10-50 ps time. Low 
excitation measurements of AR/R as a function of laser 
fluence (front and back face) show that, in this density 
regime, AR is proportional to the plasma density within 
the optical absorption depth. Hence the plasma density 
itself is decaying at this rate. 

4. In the low-fluence regime, the size of the peak of the 
reflectivity signal from the back of the film is consistent 
with a total number of electron-hole pairs in the film, 
which has not decayed substantially from the number 
originally created by optical absorption at the front face. 
Indeed, the reflectivity signal is consistent with the same 
rate of decay of the total number of electron-hole pairs as 
that observed during the 10-50 ps time range. (We note 
that this is not true in the high-flucnce regime.) 

From 1 and 2, we know that the rate of heating of the 
lattice by the plasma in the first 5 ps following excitation 
is substantially faster than 1/5 ps -1 . This is consistent 
with the electron diffraction results of Sciaini et al.— , 
which estimate that the lattice heating time, following 
photoexcitation of very thin bismuth films, is approxi- 
mately 2 ps. 

From 3 and 4 we find that, for the low-fluence regime, 
the decay rate of the total number of electron-hole pairs 
in the film is not substantially different over the first few 
picoseconds than that observed directly in the 10-50 ps 
time range. 

Our overall conclusion is that: The decay of the plasma 
density is not determined primarily by cooling of the 
plasma. Therefore, the chemical potential of the elec- 
trons and holes (near the L and T points of the Brillouin 
zone, respectively) must be different and this difference 
decays with a time-constant approximately equal to 20 
ps. The rate of equilibration of the chemical potentials 
is the electron-hole recombination rate. We find this rate 
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to be substantially smaller than that indicated by simple 
extrapolation of the low-temperature analysis of Lopez 3 
to room temperature. 

In summary, we are able to separate experimentally 
the spatiotemporal dynamics of photoexcited carriers in 
bismuth from lattice heating and strain. We measure am- 
bipolar diffusivity in the range of ^18-40 cm 2 /s and car- 
rier recombination times of ~ 12-26 ps at low excitation. 
We find evidence of rapid energy loss from the plasma 
to the lattice in both low and high-fluence excitation but 
without a corresponding decay in the plasma density af- 
ter a few ps, thus demonstrating that the electron-hole 
recombination (occurring on the order 20 ps) is not de- 
termined by the plasma cooling on this time-scale. At 
the highest density our results cannot rule out the pos- 
sibility of electron-hole recombination time of the order 
of a picosecond or less, with the plasma density being 
determined primarily by its temperature. 
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Appendix A: Electron-phonon coupling and carrier 
relaxation 

With regard to electron-hole recombination, this anal- 
ysis is similar to the approach discussed in Ref . y, Section 
V.A, but avoids the assumption of very low temperature. 
In addition, we discuss the rate of transfer of energy from 
the electronic plasma to the lattice vibration modes. It is 
important to note that the band overlap between band 5 
(valence) at T and band 6 (conduction) at L in Bi is ap- 
proximately 38 meV. Thus, the electron and hole gases at 
room temperature are in an intermediate regime between 
degenerate (Ep » kT) and classical (/(e) << 1). For 
the highly optically excited electron-hole plasma, with 
temperatures of the order of kT = 0.5 eV, the gas can 
be treated classically. In either case, the variation of the 
density of states is very substantial over the range of the 
kT and the simplifications used by Lopez are not quan- 
titatively valid. 

We allow for the possibility that the lattice tempera- 
ture 7], which determines the level of phonon occupation, 
can be different from the plasma temperature T p , which 
determines the electronic excitation in the bands. We 
also allow for the chemical potential fj,h within the hole 
(valence) band to be different from the value /i e for the 
electron (conduction) bands. For simplicity, we assume 
that the electron and hole temperatures are equal. 



More generally, we consider the scattering of electrons 
at temperature T p between two different bands, 1 and 
2, with chemical potentials [i\ and fi2, by a phonon of 
frequency ui at temperature T;. Thus, we assume that 
electronic states of energy e in band j have a probability 
of occupation, 



i 



z (e-H 3 )/kT p _|_ I : 



for j = 1 or 2. (Al) 



Let gi (e) and 92(e) be the electronic density of states 
in bands 1 and 2, respectively. The phonon occupation 
number is assumed to be determined by the lattice tem- 
perature: 



N(huj) 



1 



chu/kTi _ ^ ' 



(A2) 



We denote the electron-phonon coupling matrix element 

by 



dH 
dQ q 



(A3) 



where H is the electronic single-particle Hamiltonian and 
Q q is the normalized phonon amplitude for wave vector 
q. We will assume that the q and k dependence of M and 
LJ q is negligible for scattering between the band extrema 
1 and 2. Then, using the Fermi Golden Rule, we obtain 
a rate of scattering from states in region 1 to region 2, 
due to phonon absorption: 

Ra,l^2 = ^|M| 2 7V(M>< (A4) 
ffi(«)S2(e + MA(e)[l - / a (e+ M] de. 
The corresponding term involving phonon emission is: 
2tt 



Re,1^2 = ^\M\ 2 [N(huj) + l} X 

n 



(A5) 



gi(e)9i{e - ~h{e- hu)\ de. 

Noting that 



/i(e)JV(M[l-/a(e+M] = ex P 



Ml ~M2 
kT v 



x (A6) 



kT p kTi 



exp 
we see that 

Ra,l->2 = ex P 

Similarly, 

R a ,2^i = exp 



[l-/i(c)](JV(M + l)/a(e + M, 



Ml - A*2 


exp 


flUJ 




kTp 


kTp 


m_ 


fJ-2 - Ml" 


exp 


Tiuj 




kTp 


kTp 





Re 



ReA- 



(A7) 



(A8) 
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Then the rate of loss of carriers from band 1 (and corre- 
sponding gain in band 2) is 



dni 
~dt 



dni 



Treating the electrons and holes as classical, non- 
degenerate Fermi gases at temperature T p , with densities 
n e and n/,, respectively, we have that 



— RaS^2 — Re, 2->l _ Ra,2-tl 
fl2 — Ml ^ w 



dt 

1 — exp 



kTr, 



kT, kT n 



-Re, 1^2 



M^ Mc 

kT n 



An„ An 

! -=2 . (A13) 



exp 



Hi — /i 2 huj Tilu 



kT n 



kTr kT n 



-l>Ra 



l.(A9) 



For many of the cases of interest to us in photo-excited 
room temperature bismuth, the term in the exponent is 
small and we get 



dni 
~dt 



Ml - M2 



kTr) 



Tilu 
kfi 



kT v 



(-Ra,l->2 + Ra,2->1) 

(R aA ^2 ~ i?a,2^l).(AlO) 



The energy of one phonon is transferred from the plasma 
to the lattice for each phonon emission and vice versa for 
each phonon absorption. Thus, the rate of heat trans- 
fer from the plasma to the lattice associated with this 
scattering process is given by 



— — TlLU [R e ,l^2 + -Re,2->1 — Ra,l 

= TlLU [i? e ,1^2 — iJo,2->-l] 
7*2-/^1 . 



Tilu 



TlLU 



Tilu 



TlLU 



exp 

exp 
M2 - 



kTp 
Mi - u 2 



2 — Rajt-yl] 
TlLU [R e ^l — Ra,l 
1 





TlLU 


m~ 


kTp 


TlLU 


Ti'jj 



kT 



Mi 



kT, 



v 



kT, 



U2 



V 

Mi 



p 

TlLU 

Wi' 

TlLU 

kfi 



kT 

TlLU 



kT ri 



1 



Ra, 1^2 
Ra,2^1 



kTp 

TlLU 

' kT 



kTp 

TlbJ 



(Ra,- 



TlUJ 

kT r) 



(Ra, 



Ra, l->2 
Ra,2^1 
R a ,2->l) 
2 + -Ra,2->l) 



(All) 



The first term in the final expression of (IA11I) is the heat 
transferred to the lattice due to electron-hole recombi- 
nation and the second term is the heat due to ordinary 
heating without recombination. We see that the average 
energy transferred to the lattice per electron-hole pair 
recombined is less than the energy of one phonon: 



AE, 



rccomb 



= TlLU 



Ra, 



— Ra,2 



(RaA 



Ra,. 



< Tilu. (A12) 



Examining the expressions for i? ai i_>.2 and i? a .2->i, it is 
clear that R a ,i^2 ~ Ra,2^i when the relative change in 
the densities of states gi and occupation factors fi is small 
for a change of energy of Tilu. In this case, AE^^^, « 

TlLU. 



(As long as we are in the non-degenerate gas regime, 
minor corrections to this result can be obtained from 
the Joyce-Dixon relation between chemical potential and 
density.) We then have that 



dAn 
dt 



2An 

TlLU 

kT, 



Ra,< 



(A14) 



Tilu 
kT n 



(Ra,v—>c Ra,c—>v^) 



and we see that the recombination rate (i.e. the decay 
rate of An when T = T p ) is simply j n — 2{R a ^ v ^ c + 

Ra,c^>v) I n. 

The plasma-lattice heat transfer rate allows us to esti- 
mate the cooling rate of the plasma, considering it as a 
classical gas of electrons and holes with a heat capacity 
per unit volume equal to 3nk. Taking the above expres- 
sion for dQ/dt and allowing for scattering between mul- 
tiple bands, including intraband scattering, we can find 
the rate of cooling of the hot plasma when it transfers 
heat to the lattice: 



dTp 
dt 



TlLU 

3nfc 



E 

i,3 
TlLU 

3nfc 



Uj - /ii 



kTp 

hcu hcu 
kTi kT % 



2> 



1,3 



TlLU X 

3fc ^ 



Uj - u t 



kTn 



%0 

(T P - T) 



Ra. 



TlLU 
<$kTr) 



TlLU 

kTi 



y ^i!z^A15) 

Z ^ t, 



Thus, the plasma cooling rate (or carrier energy relax- 
ation rate) is 



7Q 



TlLU 

3>kTp 

TlLU 



Tilu 



<$kTri 



kT 

TlLU 



kT 



E 

i,3 

Et. 



Ra 



i,3 



TlLU 



<$kTrt 



TlLU 

kf~i 



(A16) 



7 P , 



where 7 P is the average carrier scattering rate, corre- 
sponding to a carrier momentum relaxation rate. The 
average phonon energy in bismuth is of the order of 10 
meV so that, for the initial, hot photo-excited plasma, 
where 3kT p — ?iw p hoton = 1.5 eV and kT = 25 meV, 
we find that the energy relaxation or plasma cooling rate 
7Q ~ 7p/375. For the final, cold (room temperature) 
plasma, T p — T and 

/ \ -I 2 



7Q 



kT 



~ 7 P 
7 P «y 



(A17) 
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It is clear that the momentum relaxation rate determines 
the decay of electronic excitations of symmetry lower 
than the full crystal symmetry. Thus, the observed re- 
laxation rate — of the order of several femtosecond — for 
the decay of the E g force in Li et al.— is consistent with a 
plasma cooling time (and associated lattice heating time) 
of a few picoseconds. 

We note that, whereas the contribution to the electron- 
hole recombination rate arises only in the energy region 
where the conduction and valence bands overlap (plus 
approximately 10 meV for the energy of the phonon par- 
ticipating in the recombination process) , the contribution 
to plasma-lattice heat transfer when the plasma is hot- 
ter than the lattice arises from all bands. In particular, 



the optical phonons throughout the Brillouin zone will 
be effective in this heat transfer. For this heat trans- 
fer, we must consider intraband scattering, where band 
1 is the same as band 2, as well as interband scattering. 
Thus, carrier recombination is strongly suppressed by the 
bottleneck in the density of states in the overlap region 
between the conduction and valence bands. In the case 
of the initially hot plasma (with typical kT p sa 0.5 eV), 
the electrons and holes are spread in energy over a range 
of energies of the order of 500 meV, whereas the overlap 
between the conduction and valence bands is only ap- 
proximately 40 meV. Moreover, the density of electronic 
states to scatter to is much smaller near the band edges 
than it is 500 meV into the band. 
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